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Abstract 

The ordinary Poisson brackets in field theory do not fulfil the Jacobi identity if bound- 
ary values are not reasonably fixed by special boundary conditions. We show that these 
brackets can be modified by adding some surface terms to lift this restriction. The new 
brackets generalize a canonical bracket considered by Lewis, Marsden, Montgomery 
and Ratiu for the free boundary problem in hydrodynamics. Our definition of Poisson 
brackets permits to treat boundary values of a field on equal footing with its internal 
values and directly estimate the brackets between both surface and volume integrals. 
This construction is applied to any local form of Poisson brackets. A prescription 
for 5-function on closed domains and a definition of the full variational derivative are 
proposed. 
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1 Introduction 



Field theory canonical formalism has some specific features which are absent in mechanics. 
The need to deal with quantities integrated over some region of space and to integrate by 
parts leeds in some cases to the appearance of surface integrals in the hamiltonian and(or) in 
Poisson brackets. Mathematicians usually prefer to consider all these surface integrals to be 
zero and construct their refined formal variational calculus by identifying integrands which 
are different by divergencies [|T| . But in field theory in some cases these surface terms are not 
zero and bear physical meaning. So, this paper is devoted to an extension of the hamiltonian 
formalism onto these divergent terms or informalizing the formal variational calculus. 

It seems quite natural to require that for closed systems boundary values should have 
equal rights with the internal ones and be determined exclusively by initial conditions to- 
gether with dynamical equations. Here we intend to decline any nondynamical boundary 
conditions, at least at the first stage of investigation. Afterwards they can be restored as 
any other constraints which can be put on the initial values of dynamical variables. We 
believe that such an approach can help us in the following papers of this series to solve some 
physical problems that are intractable by other methods. So, though the problem treated in 
this paper is a mathematical one, our motivations are physical. 

Let us first remind some important previous results on the problem. 

In the famous work by Regge and Teitelboim it was shown that to make the hamilto- 
nian dynamics of the gravitational field in asymptotically flat space existent and nontrivial 
it was necessary to include into the hamiltonian the surface integrals of special kind. This 
paper also contains, but not in a very explicit form, the acknowledgement of the physical 
meaning of surface integrals arising in the evaluation of Poisson brackets, because it was 
found that they were in correspondence with the surface terms in the hamiltonian through 
the algebra 



where 



I 5 



and popular boundary conditions are adopted for canonical variables 'jij,'^''^-^ and functions 
N (x) , M (x) , N'^ (x) , (x) . In our paper [^] it was shown that this correspondence could 
be exploited under more general boundary conditions for explicit evaluation of the surface 
terms in hamiltonian. The essence of the method was the independence of formally defined 
canonical Poisson brackets 



{F,G} 



SF SG SO SF 



6%j{x) 57r*J(x) 6'yij{x) S'ir'^^{x) 



(fx 



on the surface integrals in the hamiltonians F and G. This property of standard Poisson 
brackets is in full analogy with the one mentioned in Arnold's book on classical mechanics 
1^, where the two functions, defined only up to constants, give their Poisson bracket exactly, 
not up to a constant. So, if in mechanics the kernel of Poisson bracket consists of constants, 
in field theory the kernel of the ordinary bracket includes surface terms also. 
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In the papers devoted to the study of the Korteweg-de Vries(KdV) equation, such 
nonstandard features as the noncommutativity of variational derivatives and violation of the 
Jacobi identity were observed. In connection with this difficulties different modifications of 
the Gardner bracket were proposed 0,0, [0]- Their comparison can be found in Ref. P|. 
Unfortunately, these papers were unknown to us during our work on Ashtekar's variables 
0, where similar observations were made. Seemingly, such observations were made by 
mathematicians long before [^,[]TT|. In our paper it was conjectured that the general 
criterion for the choice of boundary conditions in the hamiltonian approach to field theory 
should be the fulfilment of the Jacobi identity for the standard Poisson bracket. 

As we recognized from the very readable book by Olver [|12[, when studying surface waves 
of the ideal fiuid Lewis, Marsden, Montgomery and Ratiu (LMMR) [|l^ proposed a modified 
form of canonical Poisson bracket^ 



{F,G} 



Sq{x) dp{x) 6q{x) 5p{x) 



+ 



an 



an 



Sq{x) 



Sq{x) 
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dS 



dS, 



where the two components of variational derivative were defined through the formula: 



r 6^F r d'' 

DgF{q,p) ■Sq= / ^ ■ Sqd^'x + f ^ ■ 5q\andS, 
Jn da Jan on 



6^F 



In 6q 



an 6q 



where partial Frechet derivative was used in the l.h.s., and through the analogous formula 
for the momentum derivative. This bracket was accompanied by boundary conditions that 
LMMR considered as necessary 



S^FS'^G S'^GS^F 



Sq 5p Sq Sp 



0. 



(1-2) 



LMMR mentioned, that earlier the hamiltonian structure for surface waves of the ideal fiuid 
for the potential fiow was discovered by Zakharov 



Here we will generalize the LMMR formula in such a way that the new bracket will fulfil 
the Jacobi identity without any boundary conditions. It will permit us to consider on formally 
equal grounds both volume and surface hamiltonians. The boundary values of hamiltonian 
variables are now obeying their own hamiltonian equations and fixation of some boundary 
conditions is simply a new constraint, that should be examined, according to the Dirac 
procedure, for the presence of secondary and higher constraints. Such a generalization of the 
LMMR bracket seems necessary also because this bracket of two differentiable functionals 
can be a functional not differentiable in the sense of ( |1 . 1|) and (|1.2| ). Therefore, in general, 
it is not possible even to check the Jacobi identity for the LMMR bracket. 



^ It is interesting to note that the LMMR work was reported at the same conference and pubhshed in the 
same journal issue as the work by Buslaev, Faddeev and Takhtajan in which the Gardner KdV bracket 
had been modified. 
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We announce more general new formulae for Poisson brackets and prove for some im- 
portant cases that they fulfil the new definition of Poisson brackets without discarding any 
surface integrals. The cases where the proofs are demonstrated are l)ultralocal bracket 
with constant structure matrix (the canonical case); 2)ultralocal bracket, depending on field 
variables but not on their derivatives (Lie-Poisson brackets are the most popular examples); 
3)nonultralocal brackets with constant structure matrix (Gardner-Zakharov-Faddeev bracket 
for KdV equation may serve as an example). 

Plan of this paper is as follows. 

In Sec. 2 we introduce necessary notations and give briefly the mathematical background: 
definitions, lemmas and formulae. Sec. 3 contains the primary mathematical motivation for 
the construction of the new brackets in the ultralocal case: the idea that Poisson bracket 
should generate the full variation of a local functional. In Sec. 4 the method for constructing 
general local brackets is presented. It is based on integration by parts of the local formula 
and gives some new proposals about handling distributions in this case. Of course, these 
calculations can be justified only "a posteriori" in Sec. 6. Sec. 5 is devoted to a definition of 
the full variational derivative as a distribution. Here we also present the most unexpected 
result of the paper: a multiplication rule for derivatives of the characteristic function. This 
rule permits to write the new Poisson brackets in the form of the old ones but with the new 
variational derivatives. Sec. 6 contains three different proofs of the Jacobi identity for the 
new brackets. The simplest proof is applicable only for ultralocal brackets with constant 
structure matrix. The more general proof is applicable to ultralocal brackets dependent on 
field variables. This proof heavily relays upon Aldersley's results [0 for higher Eulerian 
operators. The proof of Jacobi identity for nonultralocal brackets with constant coefficients 
is demonstrated too. We are sure that general proof for arbitrary local brackets along the 
lines of these partial proofs can also be constructed but it is clear, that it should be rather 
long. Probably, some other ways for this proof would be shorter. In the Conclusion we give a 
rather short resume, because a detailed comparison of our approach with Regge-Teitelboim's 
and LMMR's treatments of surface terms in the Hamiltonian formalism is postponed to next 
papers of this series. Appendix contains a collection of different ways of presentation of the 
new brackets. 

We plan to discuss the related physical problems of string theory, gauge and gravitation 
fields in further papers of this series. 



2 Notations and mathematical background 

In this paper we use the local coordinate language and instead of the manifold with a bound- 
ary consider a domain Q in having a smooth boundary dQ. The characteristic function 
of this domain is 9^ = 9{Pn), where equation Pn{x^, ...,a;") = defines the boundary. We 
do not expect that global formulation can meet with serious difficulties. 

Definition 2.1 An integral over a compact domain Q of a function of field variables 
(j)^{x),A = 1, ...jP and their partial derivatives Djcp^ up to some finite order 

F = [ d''xf{Mx),DjMx)) 
Jn 
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is called a local functional. 

All the functions / and (pA as well as their variations throughout the paper are supposed 
infinitely smooth, i.e. C°°(i?"'). We use the multi-index notations J = (ji, j„) 

Dj = ^ . ^ ^ . „ , \J\ = ji + ... + jn- 
Binomial coefficients for multi-indices are 

where ordinary binomial coefficients are 

(A = hi/ikKj-k)\) iio<k<r, 

\k J 1 otherwise. 

As the number of sums in some formulae of this paper is considerably more than ten, we 
write only one sign of summing without displaying the indices of summation. According to 
this half-Einstein rule, sum over all repeated indices should be understood. Only in those 
cases, where it is not so, we display the summation indices. Also, we do not show the 
limits of summation, because they are always natural, i.e. outside them the summand is 
simply zero. Such a nice property of binomial coefficients considerably helps us in many 
cases of changing of the orders of summation. There is also a temptation to remove the 
useless d'^x in the integrals and to write the arguments only when they can be mixed. In 
principle all integrals over finite domains would be better written over all with the help 
of characteristic function, but we will use in parallel the notions 

l^f or J Oaf or J eiPn)f. 

We denote as A the space of local functionals . It is very important that this space includes 



functionals with integrands depending on derivatives of arbitrary order [T^. Otherwise the 
Poisson brackets could go out of A. 

Definition 2.2 A bilinear operation {■, ■} such that for any F,G,H & A 

1) {F,G}eA; 

2) {F, G} = -{G, F} mod (Div); 

3) {{F, G}, H} + {{H, F}, G} + {{G, H}, F} = mod {Div); 
is called the standard field theory Poisson bracket. 

Definition 2.3 A bilinear operation {■, ■} such that for any F,G,H E A 

1) {F,G}eA; 

2) {F,G} = -{G,F}; 

3) {{F, G}, H} + {{H, F}, G} + {{G, H}, F} = 0; 
is called the new field theory Poisson bracket. 

Definition 2.4|jl2|, Definition 5.70] Higher eulerian operators Ej<^ are defined through the 
formula of full variation of local functional 

5F = Y.j^Dj[EiU)5<P^. (2.1) 
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Lemma 2.5[0, Statement 5.72] Higher eulerian operators can he given by the formula 

Of 



i^i(/) = E(-i)'^'^'''(^)^x- 



(2.2) 



Usual variational derivative (or Euler-Lagrange derivative ) is the eulerian operator of 
zeroth order. Let us mention that if J is not contained in K, then all quantities having 
multi-index {K — J) are zero. The sums over J in Eq. (|2.1|) and over K in Eq.( p.2|) are really 
finite because local functional can depend only on a finite number of derivatives according 
to Definition 2.1. 

Lemma 2.6[|1^, Statement 5.76] The eulerian operators have a property 

Ei{Djf)=Ei~\f). 



Just this property was the reason for the first appearance of these operators in paper 



IS 



Lemma 2.7||rT], Proposition 3.1] Eulerian operator of a product of two local functional 

EfUa) = E(-l)""'''"(^) {E^AU)DL-K9 + E^^{g)D,^Kf). 
Lemma 2.8(0, Theorem 2.1] Product of eulerian operators is 



^A^i(/) = E(-i)""f 

K \ 



J : 



E\'^- 



df 



Lemma 2.9[pT|, Proposition 1.1] 

Of 



Y.\ jlDK-jE-iU). 



K 



Lemma 2.10|TT|, Lemma 2.2] 

d 



df 



Let us mention that notations in [|ri[] are different because multi-indices are not used 
there and because the definition of eulerian operator Ej^ differs by factor (—1)'^'. 
We also need combinatorial identities 



Lemma 2.11|jTT|, Lemma 1.1] 

E(-i)' 

l=k 

that can be written as 

' (-m\ 



E 

l=k 



-1? 



i\k\ 



(^■_,)!(^-_A;)!(z + A;-j)! 
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Lemma 2.12|jT^, p. 616] 



^ \ / / \ ?■ - / 

1=0 \v \J 



that can be also written as 



E 



- - k)\{i - ly. t\j\{j - k)\{t - ky: 



Definition 2.13[12, Definition 5.28] The partial Frechet derivative of a function f is a 
differential operator Df^, defined for arbitrary qa as 

d 



e=0 



In our case 



The Leibnitz rule is 



I 0(Pa 



DjU9) = T.{^j}jDKfDj.K9. 



(2.3) 
(2.4) 



3 Motivation for the new brackets from the full varia- 
tion formula 

As a rule, the Poisson brackets are given by the formula 

where variational derivative is believed to be the zeroth order eulerian operator ( Euler- 
Lagrange derivative) 

and where we do not care for any surface integrals because all of them are supposed to be 
zero. Here we limit our attention to ultralocal Poisson brackets. The more general case will 
be treated in the next Section. 

Definition 3.1 The standard Poisson bracket is called ultralocal if 

{0A(a;), 0b(2/)} = lAB^ix - y), 

where the so called implectic JTSf operator Iab (or structure matrix) can depend on the field 
variables (pA^x) and their derivatives Dk(Pa{x). 



6 



These Poisson brackets, together with the local functional H, called hamiltonian, generate 
a variation of any local functional F under fixed boundary values of (pA, DjipA according to 
the formula 

6hF = {F,H} = J2 t— 
Jn 0(pA 

where 

5H(t)A = y^lABTi-- (3-1) 

The new Poisson brackets we are searching for should analogously generate for a given 
hamiltonian H a full variation of a local functional F in accordance with Eq.( p^.l| ) 



5hF = {F, H} = Y.j^ Dj (Ei{f)5H<P^ . 



where variations 6h<Pa of field variables are linearly dependent not only on E%{h), but also 
on higher eulerian operators ( p.2|) 



5H(t^A = Y.lfBE^{h). 

Evidently, jjf^ = Iab, and other coefficients i'^^b "W'ill be found below. Really, they are 
distributions, and this aspect will be treated in the next two Sections. Here we need them 
"in weak sense", i.e. as functionals defined on the standard smooth functions. We will 
show, that these coefficients can be found from the requirement of antisymmetry of Poisson 
brackets, i.e., 

E j^Dj{Ei{f)lf^EUh)) = - E j^Dj{Ei{h)lf^EUf) 

Let us consider this condition perturbatively in the order of eulerian operators \ J\ + \K\. In 
the zeroth order the antisymmetry is fulfilled due to the related property of the standard 
bracket 

r{0) _ AO) 
^AB — ^BA- 

In the first order we should have 



E E /^^(/uS^fW + E E lDj(Ei{f)if,Ei{h) 

A,B\K\=l^ A,B\J\=1-'^ ^ ^ 

= -E E / m)lf^E^{f)-Y. E / Dj{Ei{h)I%El{f)). 

A,B\K\=1^ A,B\J\=l'^^ ^ ' 

If we regroup the terms and exploit the zeroth order antisymmetry, then after relabeling 
some indices (A 5), (J ^ the above relation can be written as 

E E 

A,B \K\=V 



A,B \K\=1 ^ ^ 

EE/ DK(E'A{MlE^{h) + E'Ah)I^^lE^{f) 



A,B\K\=1' 
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Taking into account the linear independence of eulerian operators we conclude that 

EE/ ^^(/)/S^fW = E E / DjE'AMlE^ih)). (3.2) 



A,B\K\ = 1'''' A,B\K\ = 1' 



Therefore, we succeed in determining the coefficients for the \K\ = 1 case. 
Then, let us consider the next order 



EE/ Dj(Eiif)I^^lEUh)) + j: EE/ Dj(Ei{f)I^/jE^ih) 

A,B\J\=2^ ^ ^ A,B\J\=\\K\=\^ ^ ' 

+ [ E\{f)lf^EUh) = -Y. E / Dj(Ei{h)I%El{f) 

A,B\K\=2^ A,B\J\=2^ ^ ' 

-EE E [ Dj(Ei{h)lf^EUf))-Y: E / E\{h)lf^EUf). (3.3) 

A,B\J\=l\K\=l''^ ^ ^ A,B\K\=2''^ 

If we take into account the result obtained before ( |3.2| ), then the second terms in the l.h.s. 
and in the r.h.s. of equation (|3.3|) are mutually cancelled. Making the same procedure as 
was used for the first order we find 

E / E\U)lf^E^{h)= E / DJE\U)lfBEUh)). 

\K\=2^'^ \K\=2''^ ^ ^ 

So, it is clear that from the only requirement of antisymmetry we, step by step, become 
convinced that the Poisson bracket should be written as 

{F^H} = Y: j^Dj^K{EiU)lABE^{h)y (3.4) 

Now we are able to formulate 

Theorem 3.2 Formula (p.4|) gives a new Poisson bracket, if its zeroth order (\J\ = = 
\K\) term is a standard ultralocal Poisson bracket and the structure coefficients do not depend 
on the derivatives of the field variables. 

Proof. The antisymmetry is clear from the construction. Evidently, the bracket is a local 
functional and all that we are to prove is the Jacobi identity. This proof will be given in 
Sec. 6, but first there we will give a considerably more simple proof for the case when Iab 
are constants. 

Remark. It is not difficult to include the case, when I^b also depends on field derivatives, 
but it makes the proof of Jacobi identity even longer and the conditions for standard brackets 
are not so transparent. 

4 Surface terms and distributions 



The standard field theory Poisson bracket [12 



{e,g} = j:[ I E'Mi^))EU9iy)){M^),My)} 
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is a special case, which is true only under assumption that all surface terms arising when 
integrating by parts are zero, of the formula 

nJn 50(f) (x) d<j>^^\y) 



or 



where Frechet derivatives (|2.3|) are used. 

Definition 4.1 The standard field theory Poisson bracket is called local if 

{M^lMv)} = \i:[i'iB{^)Di^-iBA{y)D^t)K^-y\ (4.2) 

where the sum is of finite range in \L\. 

Usually, derivatives over only one of the two arguments are present in the formulae like 
([4 .21), because of the widely used relations 

[of -{-l)\'\D'fy{x-y) = Q, (4.3) 
I'ab = (-l)l^l+^/f 



also accompanied by 



But if not all surface terms, arising in integration by parts, are zero, then (|4.3| ) are not true. 
This observation was made by the author in |^ where it had been realized that the problem 
could be reduced to the definition of integrals like 

/ / f{x)g{y)DfD''^^5{x-y), (4.4) 

for finite domain when the test functions were nonzero on its boundary. 

The theory of distributions |TP| considers them as defined on the open space domains. 
In the literature known to us related problems of defining distributions on closed domains 
were discussed in books [^,|^, but the unique answer how to define the integral in ( [4.4| ) 
is absent there. So, here we propose a Rule which is in accordance with the results of the 
previous Section. 

Rule 4.2 

/ / fix)giy)Df^D^^^6ix - y) = [ DKfDjg. (4.5) 
Jo. Jn Jn 

This Rule is different from that proposed in |^, because the rules are compatible with 
the different Poisson brackets. 

Taken together, Eqs.([4.1|) and ( [4.5| ) give us a possibility to obtain not only the previously 



found expression (|3.4|) for ultralocal brackets, but also a more general result. Let us substitute 



3) into (O 



2^JnJn d<p^i\x) dcpf\y) 
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and exploit the Leibnitz rule (|2.4|) 



'M ^AB\-^I^L+J~M^K 



D^i^ - [^^D%^iUy)DfDtK-My{^ - y)- 



Then take off one of the integrations by Rule 4.2 



L 

M^AB 



Once more using the Leibnitz rule 



Cyb ^va 



and making changes (J ^ K), {A ^ B) in the second term we obtain 



L 

AB 



Transform here the partial derivatives into eulerian operators according to Lemma 2.9 



Dn+mIab 



X i^Dp_j+K-NE^if)DQ_K+L+J-ME2{g) -{F^G) 

make a change J J + K and estimate the sum over K according to Lemma 2.12 



E 

K 



, M [j + K [n [k [m [n [ P-J~N , 



Then we get 



ln\M \N \ P~J-N 
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X 



It is not difficult to get convinced with the help of Leibnitz rule that the obtained result 
coincides with 

lDp^Q{E^if)lABDLE%) - (F ^ G)). (4.6) 

For ultralocal case I^p = 6loIab, Iab = —Iba, and, evidently, Eq.(^.4|) can be reproduced. 
For a more general case we have 

Theorem 4.3 The new Poisson brackets, corresponding to the standard local Poisson 
brackets with constant structure matrix, are given by formula ( [4.6|) . 

Proof. The antisymmetry is evident, it is also evident that ( [4. 61 ) gives local functional. 
The Jacobi identity for this case will be proved in Sec. 6. 3. 

Remark. Evidently, the construction does not restrict us to the case /^^ = const. But 
the proof of Jacobi identity becomes more difficult in the general case. 

5 A full variational derivative and a multiplication rule: 
towards informal variational calculus 

Let us present the standard variational, or Euler-Lagrange, derivative in the form 

SF 



EMOn. 



Ha{x) 

Then it gives us a full variation 
of a local functional 

F = j enf{(PA,Dj(j)A), 
only if all surface integrals in the general formula 

SF = J2j 6^Dj{EiU)5ct^A 
are zero. 

Definition 5.1 A distribution SF/S(f)A such that in general case, i.e. for arbitrary smooth 
variations 6(f) a{x), 

^F-Zj^J^^, (5.1) 

will be called the full variational derivative of a local functional F. 

Statement 5.2 The full variational derivative can be written in the form 

j-=j:i-^y'^EiU-)Djen, (5.2) 
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where 6q is a characteristic function of the domain of integration f2. 
Proof. Through integration by parts 



= Y.j^Dj{Ei{f)5<P^. 



Statement 5.3 By using Definition 5.1 the new Poisson brackets ( |4.6| ) can he written in 
the form 

if we admit the following multiplication rule: 
Rule 5.4 

DjOiPn) X DKOiPn) = Dj+KOiPn). 

Remark. Of course, this Rule has its domain of apphcabihty only inside our procedure 
of calculating the Poisson brackets. Maybe, it can also find its place in the new theory of 



generalized functions ||22 



Proof of the Statement. Let us substitute ( ^.2|) and ( [4. 21) into (|5.3|) and take off the 
derivatives from the (5-function through integration by parts 



E(-i 



AJ\ + \K\ + \L\ 



DdDj{e{Pn{x))Ei{f)lUx)]DKe{Pn{y))E^{g) 



-DLl^DKe{Pn{y))E^{g)lUy))D.APn{x))Ei{f) 6{x-y). 

Then take off one of integrations with the help of 5-function, and afterwards use the Leibnitz 
rule 

Dj+MOaDKOfiE^ {g)DL-M[ E'^{f)I^Q 



-DK+M0nDj9nEi{f)DL^M [E^ {9)Iba 
After exploiting Rule 5.4 and integrating by parts we have 

-Ei{f)DL-M{E^{9)I^BA 

and after using the Leibnitz rule 

L \ fM' 



D,,E^{g)DL-N{E'M)lAB 
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Summing over M 

M \N 



B-i)'"'ff:.V!!l=(-i)'"*..~. (5.4) 



completes the proof by giving Eq. ( [4.6|) . 



Statement 5.5 Rule 4.2 is a corollary o/Rule 5.4. 

Proof. With the help of characteristic function 6^ we can write the l.h.s. of ( [4.5|) in the 
form of integrals over infinite space i?" 



9iPn{x))f{x)9iPniy))giy)D[^^D^^^5ix - y), 
then no surface terms arise after integration by parts and we have 

I j Df{e{Pn{x))f{x)^D^i(^{e{PM)9{y)y{^-y)- 

Remove one of the two integrations with the help of the 5 -function and obtain 
then use the Leibnitz rule 

and Rule 5.4. After one more integration by parts over i?" we obtain an integral over Q 

Again exploiting the Leibnitz rule and, afterwards, calculating the sum over M, 

K\fL + M\ ( L \ 



we obtain 



After making a change N ^ N + L and calculating the sum over L according to Lemma 
2.11 we get the r.h.s. of ( [4.5| ). The proof is completed. 
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6 Proofs of Jacobi identity 
6.1 A toy proof 

Statement 6.1.1 For constant structure matrix the Poisson bracket ( p.4| ) can be written in 
the form 

{F,G} = E JjriDfjABD,,), (6.1) 
where Df^ is Frechet derivative ( |2.3| ), and 

d f dcj 

Tr{DfjABDg^) =T.^abDj—jjjDj-jjj. 

0(Pa o<Pb 

Proof. Let us use the Leibnitz rule in Eq. (|3.4| ) 

{F,G} = Y.Iab I (^^l/^DME'MlDj^J-MEiig). 
By exploiting Lemma 2.5 it can be transformed to 



dg 

A d4>B 



Then by changing the indices M — ^ M + / and the order of summation with the help of 
Lemmas 2.11 and 2.12 we are able to estimate the sum 



ALL-K- 



As a result we obtain flO). 



Statement 6.1.2 The Poisson bracket, given by formula ( p.4| ), fulfils the Jacobi identity 
when Iab = const. 

Proof. Let us transform the expression 

dh 



{{F,G},H} = J2IabIcd I Dj- 

Jn i 



with Lemma 2.10, Leibnitz rule and antisymmetry over C ^ D 



{{F.G}.//} = i:(gQwoo 
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-D 



d(p\ 'd(p\;> 



Of ~ 



After cyclic permutation we have 

{{F, G}, H} + {{H, F}, G} + {{G, H}, F} 



E 



MJ \N 
X (dj+l 



\IabIcd f F)k+n- 
Jn 



dg dh 

d4>D 



(J) 
B 



{g^h)\ + 



where dots mean analogous terms with cyclically permuted /, g, h. By changing indices 
iV ^ + J, J ^ / + M, we get 



{{F, G}, H} + {{H, F}, G] + {{G, H], = E ^_ m) ( 



J + Nj 



I IabIcd 



X / D 



J+K+N-M- 



Drr 



99 ^ dh 

-Dj^M — — - [g ^ h) 



+ ... 



So, if we simultaneously change I ^ L, J K, M ^ —N, A ^ G and B D, then the 
expression in square brackets changes its sign while the coefficient before the brackets stands 
as itself. Therefore, the sum equals zero, and the Jacobi identity is fulfilled in this simplest 
case. 



6.2 Proof for ultralocal case 

Statement 6.2.1 Ultralocal Poisson brackets, given hy formula (|3.4] ), with the coefficients, 
depending on field variables hut not on their derivatives, exactly satisfy the Jacobi identity, 
if the corresponding standard brackets satisfy it up to total divergences. 

Proof. Let us transform the expression 

{{F,G},H} = Y: J^Dj^j(^E'^^DK-,L{E^{f)IcDEU9)))lABEi{h)^ 

according to Lemmas 2.6 and 2.7. Then, taking into account that 

dIcD 



6 



MO- 



we obtain 



E(-i) 



l/l + IE-l + ILI + l^/l 



' M ' 
J-K-L. 



j^DiJEi{h)lAB 



dl 



CD 



>M0- 



X DM+K+L-i[E^if)Ei^ig) )+DM+K+L-iIcDET{E^{f)EU9) 



(6.2) 
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Let us consider the first term in square brackets. As M = 0, then the binomial coefficient is 
not zero only when I = K + L, therefore this term becomes 



E j^Dj^K+L[lAB^-^E^{f)Eh{g)Ei{h) 



(6.3) 



After cyclic permutation of F, G, H and having in mind the symmetry in J, L we see that 
this term gives no impact on the r.h.s. of Jacobi identity if 



J dIcD , J disc , r ^^DB „ 
J-AB-;^, 1- ^AD-;^j 1- ^AC-JTl — = 



(6.4) 



But just this condition is necessary [0] for fulfilment of the Jacobi identity by standard 
Poisson brackets in this case. Therefore, we can now care only for the second term in (|67 
Once more exploit Lemma 2.7, then the term of interest is equal to 



E(-i) 



_iN|/| + |X| + |L| + |JV|| 



M 



N 



T-K-L \M, 



j^DjJlABEi{h) 



xD 



M+K+L-lJ-CD 



E^E^ {f)DN.ME'6{g) + E'2Ej^{g)DN.ME^ (/) 



If we take into account antisymmetry of the coefficient before the square bracket under 
C ^ D and its symmetry under K ^ L, we obtain 



M \(N 
K -Lj\M 



) In ^^^"^ {^abDm+k+l-iIcd 



xEiih) 



E^E^{f)D^.MEh{g) - E^E^{g)DM-MEh{f) 
After cyclic permutation of F, G, H this expression can be written as 



Y^(^_l-^\I\ + \K\ + \L\ + \N\^ 



X IabDm+k+l-iIcd 
Let us exploit the Leibnitz rule and get 

Y^(^_l'^\I\ + \K\ + \L\ + \N\ 



Ei{h)DM-ME^{g)-{G^H) 



+ .... 



' M UNUI + J^ 
J-K-L [m [ P , 



X 



'I + J - P\ fl + J - P - Q\ fl + J - P - Q - B> 



Q 



R 



S 



E>pIabDq+m+k+l-iIcdE>rE^EI^ (f) 
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X 



D 



S+N-M 



EUg)Di+j,P^Q^R.sEi{h) - (G ^ H)) 



+ .... 



Transform the coefficient before the square bracket according to Lemmas 2.8 and 2.5, i.e. 
make a substitution 

DnE^E^if) = + (^^^ _]d 



K 



' U ^ 



R+T+U-N- 



After that, it is possible to simphfy the expression before the square brackets by changing 
indices, order of summation and exphcit calculation of the four sums of binomial coefficients. 

First, make the changes T^T — K,M^M — K,N^N — K and estimate the sum 
over K according to Lemma 2.11 



E(-i) 

K 



K\ 



T 



M -K 



N -K" 



'N -T\ L + N ~r 



J< \I - K - L \M - K , 



I -L 



(during its calculation a trivial shift of argument is exploited, later we will not mention these 
details). 

Then make redefinitions Q Q + 1 — M — L, S S + M — N, R —>■ R + N and calculate 
the sum over M by exploiting Lemma 2.12 



j^\l + Q-M-L) 

fj + M + L-P-Q^ 
[ R + N 



^L + N-I\fj + M + L-P-Q-R-N^ 



N-M 

\ Q J 



S + M-N 



'I + J-P^ 
. N + R J 



'I + J-P-N-R^ 
.I + Q + S-N-L. 



After a new replacement R ^ J + L — P — Q — R — S we can estimate the sum over N 



E 

N 



' I+J-P ^ 

.N + J + L-P-Q-R-Sj 



'I-N-L + Q + S-R^ 
, I+Q+S-N-L J 



' U \/N~T\_/J + L + U-P-R\/ U ^ 
.N-T [l-L J [ Q + S-T [l-L. 



And the last summation is 



E(-i) 



' u ' 

J-L, 



'I + J-P^ 
R 



7 + J^ 



-1) 



\U\ + \L\ 



'P + R^ 



As a result, the term under study becomes 



E 



'Q + S\fP + R^ 

. Q [ P J 



X DqIcdDj+l+T+U~P-Q-R~S 



L + J ^ 
P + R-Uj 



'L + J + U-P-R^ 
Q+S-T J 



' L + J ^ 
.P + R-U , 



DpIab 



DnEi(h)DsE'i(g) - (H ^ G) 



It is not difficult to see that under the simultaneous change A^C,B^D,J<-^L,R^S, 
U ^ T, P ^ Q the square bracket changes its sign whereas the coefficient before it does 
not, i.e. all the expression equals zero. With the previous study ( |6.3|) , (|6^ ) in mind the 
proof is completed. 
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6.3 Proof for nonultralocal case 



Statement 6.3.1 Nonultralocal Poisson brackets given by formula (|4.6|) , exactly satisfy the 
Jacobi identity, when I^^ = const. 

Proof. By Lemma 2.6 and changes I ^ J, A ^ B we have 

{{F,G},H} = \Y.IcDj^Di^jl^[lf^DKEi{h)-Ei{h)I^^DK 



xE 



and by using Lemma 2.7 obtain 



I-L~M / TTi-L 



E 



I~L~M 



p ' 

L-M, 



X (EP^El:,{f)Dp_,^L+M+NE'^{g) + E^-''E^'{g)Dp_j^L+MES{f) 
Then let us exploit the symmetry L ^ M and change C ^ D 



P ^ 
L~ M, 



{{F,G},i/} = iE(-l)'™'"'^"''(^ 
X j^D,^jl^(^l'^^DKEi{h)-Ei{h)I^^DK)(^I^^^ 
xDl+m+p-i+nE^ {g) - I^c^a'^ Ec{9)Dl+m+p-iE^ {f)j^. 



Make a change P P + N in the second term 



T-L-m] '''' ^ ' [i-L-Mr'''' 



x(_i)|P|+l^l+l^l+|M| j^D,^j[{l^pDKEi{h)-Ei{h)I^^DK 

X (i?:^i?^(/)DM+L+P+iv-/i?^(^?) - (F ^ G) 
Then calculate according to the Leibnitz rule 

X: {DK-QE^AEhU)DQ+M+L+P+N-lE^^{g) - (F ^ G)) 

and, analogously, -D/+j 



{{F,G},i/} = -5: 



^ \ riV _ ^_lN|iV|riV ( P + N ^ 



:^_iyp\+\i\+m+\M\ 



'I + A fl + J - 



R J\ S J Jn 



X ( DsE'^^E^{f)Dj_R_s+M+L+P+NE'D {9) - (F ^ G) ) - I^^DnE-^ih) 



X 



H j {pS+K-QEAEc{f)Dj-R-S+M+L+N+P+QE^{g) - {F ^ G)^ 



Now we are able to sum over I according to Lemma 2.11 



^ {^_l^\P\+\L\+W\ ^ 



R + S\(J + L + M\ 
R )\R + S-p]' 



'^_l-^\P\+\L\+W\+\N\ 



'R + S\f J + L + M ^ 
, R }[r + S - P - Nj 



and obtain 



1 



{{f,g},h} = -y: 



'J + L + M^ 
,R + S-P, 



^CD ^DC\ 



J + L + M ^ 
R + S-P-Nj 



'R + S^ 
R J 



X jJ^I^^DR^KEi{h){DsE^^E^{f)Dj_n_s^M+L^P^NE^{g) - {F ^ G) 



-lLDnEi{h)Y.(^^^ 



X (j^j-R~s+M+L+N+p+QE^(g)Ds+K-QE^E^(f) - (F ^ ' 



If changes R ^ R — K are made in the first term 

J + L + M \^^r .^j ( J + L + M 



jN _ jN I 

R + S-P-K) ^^\R + S - K - P - Nj 



'R + S-K^ 
, R-K , 



X J^I^sDnEi{h)(^DsE^E^{f)Dj_n-s+M+L+P+N+KE^{g) - {F ^ G)), 
and S^S — K + Q in the second 



J + L + M 



^r + s + q-p-k)^^^ ^^^\r + s + q-k -P-Nj 



J + L + M 



X 



'R + S + Q- K\fK^ 



R 



Q J Jn 



I^AE>REi{h) 
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x(^DsE^ESif)Dj_R_s+M+L+N+P+KE^{g) -{F^ G)), 



we obtain 



{{f,g},h}=-y: 



' J + L + M ^ 
^R + S-P-K, 



J + L + M \fR + S-K\^j^ 



R + S-K-P-Nj 



R-K 



'R + S-K^ 
, R-K , 



jN jK 
^CD^AB 



DC^AB 



-E 

Q 



J + L + M \fR + S-K + Q\(K\ ^ J, 

I I „ I inni T- 



.R + S - P - K + Qj 



^CD^BA 



E 

Q 



' J + L + M ' 

.R + S-P-K-N + Q, 



'R + S-K + Q\fK\^j^ 



R 



TIM T 

^DC-'-BA 



X DnEiih) (^DsE^ES(f)Dj+M+L+p+N+K-R-sE^(9) - (F ^ G) j . 
Adding the terms with cychc permutations, we group terms hke 

DsE^ES{f)I^DnEi{h)Dj+M+L+P+N+K-R-sE^{g) - {H ^ G)), 
and, according to Lemma 2.8 substitute 

DsE^^ESif) = DsU-^ri'^V]^"-^ 



^(_1)|C^I + I-P| 



T,U 



' U ^ 



'L + T^ 



\D 



Then make a change T ^ T — L 



1 



E(-i) 



u 



,L + P-Tj 



X ( DRE-i,{h)Dj+M+L+p+N+K-R-sE^ {g) - {H ^ G) 



and S^S + L + P + N + K 



^E(-i)'^'^'^'(3( 



u ^ 

L + P-Tj 



X 



^ J + L + M ^ 
M + S + L + N, 



'R + S + L + P + N\^^ ^j. 



R-K 



'J + L + M\(R + S + L + P + N^ 



-E 

Q 



\R+S + L) ' 

J + L + M ^ 
R + Q + S + L + Nj 



R-K 



Tl\ Tl 

^CD^AB 



jN jK 
^DC^AB 



'R + S + Q + P + L + N\fK\^^ 



R 



^CD^BA 
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-^\R + Q + S + Lj\ R )\Q 

X Ds+u+T+N+K jr^) [pR^B (h) Dj+M-R-sE^ {q) - {H ^ G)\ 



So, we are able to estimate sums over P 



E(-i)'"'( 



U \fR + S + L + P + N\ 

l + p-tI[ r-k I 



E(-i)'^' 



'^\\L\+\u\+\T\ (T + R + S + N\ 
' ' \ R-K-U )' 

' U \fR + S + L + P + N + Q\ 

^l + p-t)[ r ) 



and obtain 



' J + L + M \/T + R + S + N^ 
^R + S + L + n)\ R-K-U j 



tN tK 
^CD^AB 



/J + L + M\(T + R + S + N\ ^ ^ 

[r + s + l)[ R-K-U 

J + L + M \(T + R + S + Q + N\(K\ J, 

^[r + q + s + l + n)[ R-U JUJ''''''^ 



^/ J + L + M \fT + R + Q + S + N\(K\ j^ ^ 

^\\r + q + s + l)\ R-U JUJ^^^^^^ 

xDs+u+T+N+K ^ fu)r, ,(T) \ DREB{h)Dj+M-R-sE^ (g) - {H ^ G) ] . 



Summing over L 



Y \l)\L + R + S + n) ^ ' \R + S + N + TJ' 

+^ ' \LJ\L + R + S + N + QJ ^ ' \R + S + T+N + q)' 



J + M \ 

\r + s + t + q)' 
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we get 



4 ^ 



E 



+ E 

Q 



/ J + M 
\R + S + N + T 

( J + M 
\R + S + T 

J + M 



^T + R + S + N^ 
, R-K-U , 



jN jK 
^CD-'-AB 



T + R+S + N^ 
R-K-U J 



tN tK 
^DC-'-AB 



X D 



^R + Q + S + T + Nj 

J + M ^ 
^R + Q + S + Tj 



T + R + S + Q + N\/K' 
R-U [q, 



tN jK 
^CD-'-BA 



'T + R + Q + S + N\fK^ 

R-U Jig, 



tN jK 
^DC-'-BA 



S+U+T+N+K- 



DnE'MDj+M-R^sE'riig) - {H ^ G) 



(6.5) 



Let us make change of indices S 
of the square brackets 



-S — R + J + M. We can sum over Q in the third term 



E 

Q 



J + M 



J + M-S + N + Q + T, 



'J + M-S + N + Q + T^ 
R-U 



'J + M\ J + M + K + U-R^ 



R-U . 



S-N-T 



Then after interchanging R ^ S , J ^ M, B^D,A^C,N<-^K and U ^ T we see 
that the first term in the square brackets stands as itself, the new second term is equal to 
the old third and vice versa. The fourth term transforms into itselfl 



Q 



J + M \ 
J + M-S + Q + T) 

f J + M ^ 

[j + M -R + Q + U, 



'J + M-S + Q + T + N\/K\ 
R-U )[q) 

fJ + M -R + Q + U + K\fN\ 
\ S-T )\q]- 



Evidently, the round bracket in Eq. (|6.5|) changes its sign, so the expression is zero and the 
proof is completed. 



7 Conclusion 

It is clear that the above results can be applied to the field theory on manifolds with a 
boundary simply by postulating the Rule 4.2 and independently of any reasoning about 
characteristic functions. The new Poisson structure permits to consider dynamical problems 
in which boundary values of hamiltonian variables are treated on equal footing with their 
internal values. The dynamics of field variables on the boundary is determined by both 
volume and surface parts of the hamiltonian. Any choice of boundary conditions is in fact a 

•^We are able only to verify this fact by computer simulation. 
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constraint in the phase space and should be treated along with standard procedure of search- 
ing for secondary and higher constraints. These boundary conditions do not interfere with 
the dynamical equations inside the domain until we start solving elliptic type constraints, 
such as the Gauss law in gauge theories. Then the nonlocal dependence appears, including 
the dependence of surface variables, and surface part of the hamiltonian begins to influence 
the equations of internal variables ( "divergencies cease to be divergencies" in terminology of 
Arnowitt, Deser and Misner P^, p. 434] ). 
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Appendix 

Here we list different forms in which the new Poisson brackets can be written for the 
local case ( |4.2| ): 

1) through the full {not standard) variational derivatives defined by the formula ( ^.21) and 
with account for Rule 5.4: 

2) through higher eulerian operators (|2.2| ): 
where 

N 

3) through Frechet derivatives (|2.3| ): 

4) through some matrix notations : 

defined below 

r - D / 

^JK,LA^,AB ~ \-^\^J+K~L-MJ-AB- 



25 



